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Abstract 

It has been suggested that the volume dependence of the spectral weight could be 
utilized to distinguish single and multi-particle states in Monte Carlo simulations. 
In a recent study using a solvable model, the Lee model, we found that this criteria 
is applicable only for stable particles and narrow resonances, not for the broad 
resonances. In this paper, the same question is addressed within the finite size 
formalism outlined by Liischer. Using a quantum mechanical scattering model, the 
conclusion that was found in previous Lee model study is recovered. Then, following 
similar arguments as in Liischer's, it is argued that the result is valid for a general 
massive quantum field theory under the same conditions as the Liischer's formulae. 
Using the spectral weight function, a possibility of extracting resonance parameters 
is also pointed out. 
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1 Introduction 



Low-energy hadron-hadron scattering plays an important role in the under- 
standing of non-perturbative physics of strong interaction. Due to its genuine 
non-perturbative nature, such problems can only be studied from first princi- 
ples using non-perturbative methods like lattice QCD. Liischer has outlined a 
finite-size formalism which enables us to calculate the elastic scattering phase 
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shifts using lattice simulations [T|2|3Pf5] . Over the years, extensive numerical 
simulations have been carried out to the study on hadron-hadron scattering us- 
ing Liischer's formalism, both within the quenched approximations and using 
gauge field configurations with dynamical quarks [6.7,8]l^pUllllfl2fl3lll4fl5fl6j . 

In lattice study on hadron spectroscopy and hadron-hadron scattering, the 
most important physical quantity is the energy of the system which is ob- 
tained via the measurements of various correlation functions. However, since 
a quantum field theory does not conserve particle numbers in general, the 
distinction between single- and multi-particle states becomes an important 
and delicate issue in lattice calculations. In the infinite volume, the differ- 
ence is obvious since they have different kinematic behaviors: single-particle 
states have discrete energy eigenvalues when viewed in their rest frame while 
multi-particle states usually have continuous spectrum starting from the corre- 
sponding threshold. However, when performing a lattice simulation in a finite 
volume, all energy eigenvalues in the finite box become discrete. Therefore, 
other means have to to be applied in order to identify the particle nature of a 
corresponding state. 

In principle, differences between single- and multi-particle states still persist 
in a finite volume. For example, although both have discrete spectra, the level 
spacing between neighboring multi-particle scattering states becomes infinites- 
imally small while the level spacing between the neighboring single-particle 
states remains finite as the volume goes to infinity. However, it is difficult 
to utilize this difference as a practical criteria since this requires the com- 
putation of excited energy eigenvalues in Monte Carlo simulations which is 
usually quite challenging. Another method suggested by various authors is to 
use the so-called spectral weight as the identifier. This is the quantity which 
can be measured directly (and relatively easily) from Monte Carlo simula- 
tions, together with the corresponding energy eigenvalue. In a finite volume, 
the volume dependence of the spectral weight for a eigenstate is expected to 
show different behavior for single- and multi-particle states. For example, one 
expects the following empirical rule: the spectral weight to show little volume 
dependence for a single particle state (if properly normalized), while for a two- 
particle state, it is expected to show a dependence where L being the 
size of the cubic box. This expected difference in volume dependence can be 
measured in lattice simulations by performing the same calculation in two dis- 
tinct volumes. As an example, this strategy has been used in Ref. p/Tj to study 
the possible penta-quark state. Using this technique, the authors concluded 
that the expected penta-quark (single-particle) states measured in their lattice 
calculations are in fact kaon-nucleon two-particle scattering states. However, 
this conclusion is not so settled even in the first-principle lattice QCD cal- 
culations [T8|19f20j . Therefore, the volume dependence of the spectral weight 
indeed can provide us useful information about the particle nature of the cor- 
responding state. 
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In a previous model study, we have shown that the above mentioned empir- 
ical rule to distinguish single- and multi-particle states are in fact only valid 
for stable particles and narrow resonances. Using a solvable model, the Lee 
model, we showed that this rule is violated for broad resonances [2T]. A general 
formula for the spectral weight was obtained which can show either single- or 
two-particle volume behavior depending whether the width of the resonance 
is narrow or broad. 

In this paper, we attempt to generalize this conclusion that we obtained in the 
Lee model, to the case of general massive quantum field theory. For this pur- 
pose, the general Liischer's formalism is adopted. In previous studies, people 
have been focusing mainly on the energy eigenvalue (which directly enters the 
famous Liischer's formula) of the system within Liischer's formalism. However, 
since the spectral weight W{E, L) of a given state is intimately related to the 
overlap of the exact energy eigenfunction with the free scattering states, we 
have to study the wavefunction of a energy eigenstate in a finite volume. In this 
paper, our study focuses on the wavefunction in the Af sector and a formula 
for the spectral weight is thus obtained within the non-relativistic quantum 
mechanics model. By studying the volume dependence of the spectral weight 
in the large volume limit, we arrive at the same conclusion as we drew from 
the previous Lee model study. Then, following Liischer's arguments, this re- 
sult is generalized to massive quantum field theory. Our results also show a 
possibility of extracting the resonance parameters from the spectral weight 
function on various volumes. 

This paper is organized as follows. In Sec. we briefly review the quantum- 
mechanical model in the infinite volume. In Sec. [3l the quantum-mechanical 
model is studied on a three-dimensional torus of size L. In this section, we de- 
rive the relevant formulae for the spectral weight function and study its volume 
dependence. It is found that similar conclusion is reached as in our previous 
study using the Lee model. We then argue that, under the same restrictions as 
in Liischer's formula, our results found in the quantum-mechanical model can 
be generalized to massive quantum field theory. The possibility of extracting 
resonance parameters from spectral weight is also discussed. In Sec. IH we will 
conclude with some general remarks. Details on the evaluation of a function 
F(fc^) are listed in the appendix. 



2 The Model in the Infinite Volume 

Consider a quantum mechanical model whose Hamiltonian is given by: 
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where the potential V{r) is zero for r > a with some a > 0. We now dis- 
cuss the energy eigenstates satisfying: i7\l/(r) = i?\E'(r). One can expand the 
eigenfunction in terms spherical harmonics: 

^(r) =^;„(r)F,^(n) . (2) 



with: r = rn and ipimi'f') is the radial wave-function satisfying the radial 
Schrodinger equation: 

f d'^ 2 d 1(1 + 1) . ... x\ , . X 

\T^ + -:; —r^ + ^ - 2mV{r) ^;„(r) = . (3) 

\ dr^ r dr I 



where E = k'^/{2m) being the energy eigenvalue of the state. It is well-known 
that, there exist only one solution to the radial Schrodinger equation that is 
bounded near the origin. This solution will be denoted as: ui{r; k). To fix the 
normalization, we impose the condition: 

lirn r~'M;(r; k) = 1 , (4) 



and the solution to the radial Schrodinger equation then has the form: 

'4'imir) = biraUi{r; k) , (5) 

with some constant him to be fixed by other conditions (normalization, bound- 
ary conditions, etc.). 

In the region r > a where the interaction vanishes, the solution M;(r; k) are 
expanded in terms of spherical Bessel functions: 

ui{r] k) = ai{k)ji{kr) + (3i{k)ni{kr) . (6) 



The coefficients ai{k) and (3i{k) have simple relation with the scattering phase 
shift: 

tanMfc) = ^. (7) 
ai{k) - il3i{k) ' ai{k) ^ ' 



In the low-energy limit: A; ^ 0, one normally defines: 

a^i = \mik^ai{k) , A° = hm fc-'-^Al'^) , 
fc— >o fe— »o 



^ In this paper, we have adopted the same convention as in Ref. [1] which agrees 
with Messiah's book [22]. 
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and the threshold parameters: 



(9) 

a. 



In particular, ao for / = is referred to as the s-wave scattering length. Other 
a^'s for / > are sometimes also called scattering lengths in the corresponding 
channel, although they do not have the dimension of a length 0- The threshold 
parameters ai are important because they characterize the behaviors in low- 
energy scattering processes. For example, we have: 

5i{k) ~ a^pi+i ^ , (mod vr) . (10) 



3 The Model on a Torus 



We now enclose the system we discussed in the previous section in a large cubic 
box and impose the periodic boundary condition in all three spatial directions. 
The potential itself is also modified to Vl{y) by periodically extending over 
the whole space: Viijc) = J2nez'^ ^(k + ^M)- For later convenience, we define 
the the so-called "outer region" as: 

f] = {r: |r + nL|>a, for all n G Z^} . (11) 



This is the region where the potential vanishes identically. We assume the 
size of the box is L which is much larger than any of the physical scale in 
the system. In particular, we need to have L ^ 2a so that the outer region 
admits free spherical wave solutions (asymptotic states). We now would like to 
study the change in the energy eigenvalues, the corresponding wave-functions 
and their possible connections with the scattering phase shifts in the infinite 
volume. Our discussion here will focus on the case of a cubic box whose relevant 
symmetry group being the cubic group 0(Z). Generalization to an arbitrary 
rectangular box can be performed easily by changing the symmetry group to 
the corresponding ones {D4 or D2, etc.). 

Since the boundary condition breaks rotational symmetry explicitly, we an- 
ticipate that energy eigenstates of the system will not have a definite angular 
momentum in general. To be specific, the original eigenstate in the s-wave will 
acquire mixtures from higher angular momentum modes (mainly / = 4 for a 
cubic box). However, since the original radial wave-function Ui{r;k) and the 

^ Assuming a|' 7^ which is usually the case. 

Prom normalization condition (jH), it is easy to verify that the spectral parameters 
ai has the length dimension of 2/ + 1. 
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spherical harmonics forms a complete set in the functional space, we may still 
expand the true eigenfunction in the box in terms of them: 

^(r; k) = J2 hmUi{r- fc)F,™(n) . (12) 

Im 



where the coefficients are to be determined by boundary conditions and nor- 
malization. 

In the outer region f2, the solution are those singular, periodic solutions for 
the Helmholtz equation. Thus we may write: 

^{v-k)l^n = Y,vi^Gi^{v-e) . (13) 

Im 



In the meantime, the outer solution can also be expanded in terms of spherical 
harmonics and the spherical Bessel functions ji{kr) and ni{kr): 



An 



Yim{^r)ni{kr) + ^ Mlm;l'm'Yl>m'{^r)jl'{kr 



I'm' 



,(14) 



The explicit expression for A4i'm';imi^i) given in Ref. [1] which we quote 
here: 



Mim-Ak') = E ^^^^Sfe#^V(2/ + + l)(2j + 1) 



7r3/2g«' + l 



X 



I' J 

^0 



m m' —s 



(15) 



Here we have used the Wigner's 3j-symbols and q = kL /{2n). The zeta func- 
tion Zim{s,q'^) is defined as: 



yimin) 



r - q^y 



(16) 



From the analytically continued formula, it is obvious from the symmetry of 
0(Z) that, for / < 4, the only non- vanishing zeta functions at s — 1 are: Zqq, 
and Z4^Q. This is in accordance with the fact that s-wave and gf-wave mixes 
with each other in a cubic box. In what follows, we will focus on the s-wave 
eigenfunction. 
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3.1 Liischer's formula in the Af sector revisited 



In the remaining part of this paper, we will be only concerned with the energy 
eigen-functions in the Af sector, which is the analogue of s-wave in a cubic 
box. 

A good approximation for the s-wave dominated eigenfunction can be written 
as a superposition of / = and / = 4 spherical harmonics with the s-wave 
component much larger than that of g-wave. To explicitly construct this type 
of wave-functions, we notice that the eigen-function in Af sector has to be 
invariant under cubic symmetries. It is easy to verify that, there are only two 
homogeneous harmonic polynomials which are invariant under cubic symmetry 
up to / < 4. They can be conveniently expressed as: 



3^00 = ^ , 3^40 + ^(3^4,4 + 3^4,-4) = ^(x' + y' + z'- ) .(17) 
V47r 14 4v7r V 5 

So, we may write the eigen-function in Af sector as: 

^(^^)(r; k) = boouoir; k)Y,o + h^oU,{r- k) (^^o + ^(^^M + l4,_4)j 



with 1 640 1 ^ 600 in the large volume limit. In other words, to ensure cubic 
symmetry, the general coefficients him at / = 4 with different m values must 
have definite ratios. In the outer region, using relation ([6]), we have: 



^^(^^)(r; fc)|ren = &oo[«oJo(A;r) + /3ono(fcr)]Foo(a) 

+ 64oKj4(A;r) + ^^riiikr)] f^o + ^(^4,4 + ^,-4)] + {iO) 



:i^ 



On the other hand, we know that, in the outer region f2, the eigen-function 
can also be expanded into singular periodic solutions of Helmholtz equation. 
Since Gim = 3^«m(V)G(r; k"^) with G{r; k"^) being rotationally invariant, we 
see that in order to keep the eigen-function invariant under cubic symmetry, 
we must have the combination: G/^q + -\/70/14(G'4^4 + G'4,_4) in the expansion. 
Thus we may write: 

^(^'Hr;A:)Uc= (^) ^^00 

^ For simplicity of the following equations, we have scaled out an overall factor 
{4:71 /k) and an extra factor of for the coefficient of G40. 



Goo + f G40 



70 
IT 



[6*4,4 + G'4,_4] + 



.(20) 
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The fact that such a combination respects cubic symmetry can also be checked 
exphcitly. Using the expressions (fTi|) and (flSll . we may write the expansion 
for Goo as: 



00 



k 

4:71 



(no + moojo)^oo + ]j^mo4ji + "^[^44 + Yi-ijj 



(21) 



where we have introduced: moo = -^oo;00 and mo4 = 2y 3/77Vl40;oo for later 
convenience (see Ref. |1] for the notation). Similarly, for the higher angular 
momentum functions, we have: 



k^ 

G^o = — [^,41^40 + A^40;00j0^00 + •^40;20j2^20 
+ A440;40j4^40 + A^40;44 ^4 (^44 + >4 -4)] , 

k^ 

G4_4 + G4._4 = — [^.4(14^4 + ^4,-4) + 2A^44;00j0^00 + 2A^44;20j2^20 

47r 

+ 2A^44;40j4nO + (A^44;4,-4 + A<44;44) J4(n4 + ^1,-4)] , (22) 

In the above expansions, we have also utilized the following properties of the 
matrix elements Mim-^i'm'- 



M 



Imd'm' 



Mvrr,' 



I'm' dm 



M 



l,—m;l',—m' 



(23) 



Note that in the expansion of G40 and G44 + ^4^.4 in Eq. fl22|) . there are 
terms with / = 2, m = spherical harmonics. However, when we construct 
the combination G40 + (\/70/14)(G44 + (74^-4), the terms with / = 2 cancel out 
explicitly since: 7W40;2o + (\/70/7)A^44;2o = which can be checked by looking 
into Table E.l in Ref. [1]. Therefore we finally have: 



G^4o H — T-riGu + ^4 4) — — 
14 Att 




( ^70 

+ (n4 + m44j4) ( + ^^[^44 + '5^4,-4] 



(24) 



where 77144 = -^40;40 + ■ ■ ■• At this stage, it is worthwhile to point out that, 
^^00) '"^04 and 77744 that we introduced here are exactly those reduced matrix 
elements of M. in the sector. Please refer to Ref. [1] for further detailed 
explanations (especially Table E.l and Table E.2 in the reference). 

Collecting relevant information from the expansions obtained thus far, i.e. 
Eq. ([20]), Eq. ^ and Eq. ([21]), we have: 
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+ I \l ^moi + v^omu 



no + moo jo + ^"40 Jy"^04jo ) ^oo 



V 70 

J4 + W40'^4 I {Y40 + -^(^4,4 + ^4-4) 



■(•2-5) 



We should now match the two solutions given by Eq. f|T9l) and Eq. (I25l) in the 
outer region Q. This yields the following set of linear equations: 



Woo = boopo , ^'oo "^00 + 



yt;40"^04 I = &00"0 , (26) 



^^00^^40 = hoPi , Voo\ J -^mo4 + ^^40^44 I = &4oa4 ■ (27) 



These four equations can be viewed as a set of homogeneous linear equations 
for the four coefficients: foo, &005 VooV^q and 640. Demanding a non-trivial so- 
lution to exist requires the corresponding determinant of the 4x4 matrix to 
vanish. Another simple way to proceed is to divide the second equation by the 
first and similarly divide the fourth one by the third. This will eliminate all 
coefficients except for t>4o. We then arrive at: 

[12 nr 

cot (5^°) = moo + \/ ^^^4o"^04 , cot 6^'^^ = tum + J —mo^/v^o . (28) 



Eliminating f4o from the above two equations then yields: 

(cot5(°) - moo) (cot(5(^) - m44) = mo4mo4 • (29) 



This is exactly the equation obtained by general Liischer's method when we 
only consider the mixing between / = and / = 4 waves [1]. Therefore, using 
more explicit construction, not only have we recovered Liischer's formula, we 
also obtained an explicit approximate expression for the energy eigen-function 
in the Af channel which is given by Eq. f|T8l) in general and given by Eq. fl20l) 
in the outer region. 



3.2 The spectral weight function and its normalization 



Now we would like to derive a formula for the spectral weight function which 
can be measured in a Monte Carlo simulation. Instead of working with gen- 
eral states, we will focus on the single- and two-particle states. These states 
naturally arise in the lattice study of hadron-hadron scattering and hadron 
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spectrum. In such simulations, one constructs an operator (also known as the 
interpolating field operator), or operators if more than one is needed, within 
a specific symmetry sector of the theory. The correlation matrix among these 
operators are then computed by ensemble averaging over different gauge field 
configurations that is generated in a Monte Carlo simulation. 

For this purpose, we pass over to the second-quantized version of our quan- 
tum mechanical scattering model. In this model, two distinguishable particles 
scatter via a potential V{r) where r being the distance between them. The 
center-of-mass coordinate of the two-particle system is separated out and the 
mass parameter m in the Hamiltonian ([T]) refers to the reduced mass of the 
two-particle system. For each type of particle, a local scalar field operator 
7rj(x, t), with i = 1,2 designating different types of particles, is introduced 
together with its momentum space counterpart: 

7r,(x, t) = -1= ^ ^,(p, t)e^P- , n,ip, t) = ^ J rf'x7r,(x, t)e-'^-^ (30) 



They satisfy the usual equal-time commutation relations: [7rj(p, t), 7rj(k, t)] = 
^pk^ij- Using free states made up of two particles, one from each type, one can 
form a state: 

1$) = Ot(0)|0) = -^E"^(P)^I(P'0)vf|(-P,0)|0) , (31) 



with the interpolating operator 0{t) defined by: 

^it) = ^E**(P)^i(P'^)^2(-P,t) . (32) 



Requiring such a state to be normalized as: ($|$) = 1 yields the condition: 

^EI'^(P)I' = 1- (33) 
^ p 

If such a state were a bound state of two particles, $(P) would be the corre- 
sponding momentum-space wavefunction normalized according to the above 
equation. 

We can now define the corresponding correlation function: 

C{t) = (O|0(t)Ot(o)|O) = Y.mO^iOmi'e-''' , (34) 

E 



For simplicity, we assume that the two particles are distinguishable. 
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where E and \E) represents the eigenvalue and eigenstate of the full Hamilto- 
nian, respectively. By fitting the time- dependence of the correlation function 
obtained from Monte Carlo simulations, the exact eigenvalue E, and the cor- 
responding spectral weight function W{E), which is the coefficient in front of 
the exponential, is obtained. If we denote the overlap of two wavef unctions: 

0{E) = {E\O\0)\0) = J dh,dh2{E\n,r2){n,T2\O\0)\0) , (35) 

the spectral weight function is simply given by: 

W{E) = |(E|Ct(0)|0)p = \0{E)f . (36) 

At this point, it is worthwhile to point out that the spectral weight function 
W{E) defined above depends explicitly on the normalization of O. 

Due to translational symmetry, the exact wave-function (ri,r2|-E) will only 
depend on the relative coordinate r = r2 — ri. It is independent of the center- 
of-mass coordinate r^. This means that, if the eigenstate \E) is normalized 
according to {E\E) = 1 as it should, the wave-function (ri,r2|-E') = 
should be normalized according to: 

Jdh\{v\E)\^ = ld'r\¥''t\v;k)\' = ^. (37) 

Therefore, in order to compute the volume dependence of the spectral weight 
function, we first have to fix the normalization of ^'^'^i •*(r; k) according to this 
convention. 



3.3 Normalization of the energy eigenstates in sector 

As discussed in the previous subsection, the wavefunction in the Af sector in 
Eq. (|T8|) must be normalized properly on the torus % according to Eq. fl37|) . 
The integral of the eigen-function on the torus runs over two regions: the inner 
region where the explicit form of the wavefunction is not known and the outer 
region Q where an approximate form of the function is given by Eq. ( l20l) . 
Although we do not know the exact form of the eigen-function in the inner 
region, we do know that the eigenfunction is bounded in this region. Since it 
is assumed that the interaction region is of size a with a <^ L, therefore the 
integral in the normalization condition (1371) is dominated by the integral of the 
function in the outer region Q. Therefore, we may modify the normalization 
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condition to: 



Since in the large volume limit, the eigen-function is dominated by the s-wave 
contribution, we may use the first term in Eq. ( l20l) and write: 



2 



l^^ool 



I dh\Gooir;k)\'- J dh\Gooir; k)\^ 

\T3 B J 



7? ' ^''^ 



where the second integral is over the interaction ball region: B = {r : r < 
a, mod L}. We now use the definition for Goo: 



where the summation of p = (27r/L)n is for all three-dimensional integers: 
n G Z'^. Substituting this expression into the first term and Eq. fl25l) into the 
second integral in Eq. ( l39l) we get: 

k"" 1^1 k^ " 



E / 2 _ fe2^2 - J ^'^'^ (^o(A:r) + mooJo(A;r)) .(41) 



16tt^\voq\''L^ AttL^ ^ (p2 - A;2)2 IGtt^ 



The integral in the second term maybe evaluated directly within r < a. We 
thus obtain: 



1 4.^ 1 



\vqq\^L^ k^L^ ^ (p2 - k^f 2A;2 sin^ A 



/ sin ka\ 

1- — h°^<*''+2^' 



.(42) 



where we have utilized the definition: moo = cot A. In the large volume limit, 
the first term on the r.h.s. of the above equation is much larger than the second 
(see appendix |A] for the explanation of this assertion). If we drop the second 
term, we then arrive at: 

47r\2 , ,„ „ / 1 ^ 1 47r 

^) KPL»:.4.^^i:^-j-^j (43) 

where we have defined the function: 
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where we have introduced a cutoff function /(p^) to regulate possible ultra- 
violet divergences. The property of this function in the large volume limit is 
addressed in appendix |X1 The relevant formula for us is given by Eq. (lA.lOp . 



3.4 Spectral weight in Af sector 



We now evaluate the spectral weight using Eq. fl36|l with the exact energy 
eigen-function given approximately by: \E'*-'^i ^(r; k) ~ (47r/A;)fooG'oo(r; k). The 
overlap of the two wave-function is approximately given by: 



Using the expression ( 143|) and the expression in Eq. (lA.lOj) . we finally obtain 

W{E) = |0|2 as: 

W(E) = 8^fc|v^^(fe')P ^ 87rk\Mk')\' .^g. 

^ ' cot6o{k) + ^csc^So{k) cot So{k) + ^ csc^ 6o{k) ' ^ ^ 



where the function ipL{k'^) is defined as: 



In the large volume limit, following similar derivation as in our discussion of 
function F{k'^), this function goes over to: 



Thus the function ipilk"^) has little volume dependence in the large volume 
limit. Therefore, the explicit volume dependence of the spectral weight func- 
tion W{E) comes mainly from the denominator in Eq. (H6!) . Normally, if 
cot5o(A;) is not changing rapidly, the second term in the denominator of 
Eq. (H6l) . which is proportional to L^, dominates the result and one finds that 
the spectral weight is proportional to This is the typical two-particle 

spectral weight function. However, if there exists a rather narrow resonance at 
energy E = E^, then close to this resonance energy, one has approximately: 

cot5(E)^:^^, (49) 
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where T is the physical width of the resonance. In this case, we obtain: 



W{E) 



(50) 



{E, -E) + ttE, 



1 r 

'*AE 



If r/AE -C 1, then the quantity in the denominator is dominated by the first 
term and the spectral weight shows a typical single-particle behavior. This 
means that an extremely narrow resonance behaves like a stable particle. If 
on the other hand T/AE ^ 1, which is always true for an extremely large 
volume (assuming the width of the resonance remains finite), the denominator 
is dominated by the second term and the spectral weight itself is roughly 
proportional to which is typical for a two-particle scattering state. We 
therefore arrive at the conclusion that the volume dependence of the spectral 
weight near a resonance is controlled by the ratio (T/AE). 



3. 5 Generalization to massive quantum field theory 

Our results on the volume dependence of the spectral weight is obtained within 
a quantum mechanical model. In this subsection, we would like to generalize 
these results to massive quantum field theory, following the line of arguments in 
Liischer's formalism [5] . Using an effective Schrodinger equation (derived from 
the Bethe-Salpeper equation) [2], Liischer has argued that, if the size of the 
box is large enough such that all quantum field theory effects are suppressed 
exponentially, the results obtained within the quantum-mechanical model can 
be carried over to the case of massive quantum field theory literally [2115] . Here, 
we will assume that the same conditions are satisfied and thus our results 
obtained within the quantum-mechanical model are expected to be valid for 
massive quantum field theory. 



3. 6 Possibility of extracting the resonance parameters from the spectral weight 

The relation established in Eq. (1501) opens up a possibility for extracting the 
width of a resonance if the spectral weight can be measured in the simulation. 
Assuming that there exists a single resonance in the energy region that we are 
interested in, and the contribution from this single resonance dominates the 
scattering, we simply rewrite Eq. fl50|) as: 




1 1 




(51) 
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Therefore, by fitting the function 1/W{E, L) for different E and L (hence 
different AE as well), it is possible to extract the width parameter F together 
with the resonance position E^ of the resonance. Note that in previous lattice 
calculations, focus has been mainly put on the energy levels, i.e. the values 
of E, only. No attention is paid to the associated spectral weight function 
W{E, L) which in fact can be obtained from the fitting procedure of the cor- 
responding correlation functions with almost no extra costs. The study in this 
paper indicates that, the spectral weight function at various volumes also con- 
tains valuable information about the scattering and might also be utilized in 
some way. In fact, it can be used as an cross-check for the scattering phase 
obtained from the energy levels. Of course, this is only a possibility at this 
stage. The feasibility of this method has to be check in realistic simulations. 



4 Conclusions 

In this paper, we have studied the volume dependence of the spectral weight 
function which is accessible in Monte Carlo lattice simulations. Motivated by 
our previous study in the Lee model, it is expected that the spectral weight 
function shows little volume dependence for a stable or narrow resonance while 
for a broad resonance, it exhibits a typical dependence, the same as a 
two-particle scattering state. To verify this scenario, Liischer's formalism is 
adopted. It is first shown in a quantum mechanical model and then generalized 
to any massive quantum field theory, assuming that the polarization effects 
are exponentially suppressed following Liischer's arguments. In particular, we 
expect this scenario to be true also for QCD which governs the scattering of 
hadrons and therefore our result is relevant for lattice QCD simulations. Our 
final result for the spectral weight is summarized in Eq. (I46p which exhibits 
either single- or two-particle volume dependence depending on the value of 
r/AE where F is the physical width of the resonance and AE is the typical 
level spacing near the resonance in the finite volume. Possibilities of using the 
formula to extract the width of a resonance is also discussed. 
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A The function F{k'^) 



To study the normahzation of the wavefunction ^("^i )(r; k) in the large volume 
limit, we define the function: 



where we have introduced a cutoff function /(p^)- The relevant function ap- 
pearing in the normalization condition (l43|l is given by the derivative of F{k'^) 
with respect to k"^: 



We now follow the argument in Ref. [21] to estimate the value of F{k'^) for 
arbitrary value of fc^ in the large L limit. We separate the summation into two 
parts with: |p^ — > e and |p^ — < e. The first part goes smoothly to the 
principle-valued integral while the second summation may be written 

as: 



- y 



p,|p2-/c2|<e 



P^ 



A;2 L3 



- y 

T3 ^ 



1 



n=— oo •f* 



p^ + nAp2 — 



TT 



cot 



TT 



. AP^ . 



(A.3) 



where p^ is the value of p^ that is closest to A;^; Ap^ is the typical level spacing 
between neighboring p^ values which can be estimated by: 



(M 



r27r Jp2Ap2 = 1 ^^ L^Ap^ 



(270^ 



(A.4) 



Therefore we obtain: 



cot 



TC 



, Ap2 



(A.5) 



However, since it is easy to verify that: 

goo(i;g^) ^ k 



i^(fc^) = %^-^cot5o(A:), (A.6) 
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where we have utihzed the approximate relation (Liischer's formula): 

2oo(i;g') 



cot 6o{k) 



(A.7) 



We therefore seem to have: (pik"^) = in which case we recover the DeWitt's 
formula: 



6o{k) 



-TX 



(A.8) 



If one evaluate 0(/c^) explicitly, one gets: 
c/3p 1 



ct>{k' 



V 



(27r)3 p2 - P 



47rA + 277 /c In 



k-k 



K + k 



(A.9) 



with a sharp momentum cutoff A. This expression indeed goes to zero if we 
drop the constant term and taking A — oo. Consequently we have for the 
function F'ik"^): 



cot 



Ap2 



+ 



%Txk [ \ Ap2 )\ 4Ap^ 

1 k 

cot 5o(^) H 7 — ^ csc^ ^oik) 



CSC 



TT 



. Ap2 



4Ap2 



(A. 10) 



where in the second line we have used DeWitt's formula. Since Ap^ oc L^^, 
we find that F'{k'^) cx in the large volume limit. This justifies the assertion 
made after Eq. fH2|) in the main text. 
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